AD-A1 12  122  ARMY  ELECTRONICS  RESEARCH  AND  DEVELOPMENT  COMMAND  FO— ETC  F/S  20/5 
QUANTUM  MECHANICS  APPLIED  TO  DYE  LASERS. (U> 

SEP  81  C  K  PEARCE 

UNCLASSIFIED  DELNV-TR-0024  NL 


ADA  11 2122 


A*#*  v . 


Report  DELNV-TR-0024 


QUANTUM  MECHANICS  APPLIED  TO  DYE  LASERS 


by 

Carol  K.  Pearce 


September  1981 


Approved  for  public  release;  distribution  unlimited. 


MAR  i  8  1982 


U.S.  ARMY  ELECTRONICS  RAO  COMMAND 

NIGHT  VISION  &  ELECTRO-OPTICS  LABORATORY 

FT.  RELVOIR.  VIRGINIA  22000 


Destroy  this  report  when  it  is  no  longer  needed. 
Do  not  return  it  to  the  originator. 


The  citetion  in  this  report  of  trade  names  of 
commercially  available  products  does  not  constitute 
official  endorsement  or  approval  of  the  use  of  such 
products. 


i 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  of  THIS  PAGE  (When  Date  Entered) 


REPORT  DOCUMENTATION  PAGE 


■  REPORT  NUMBER 

1DBLNV-TR-0024 

♦.  title  (end  Subtitle) 


ft 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


CIPIENT'S  catalog  number 


3ECIPIENT*S  CA 


s.  TYPE  OF  REPORT  &  PERIOD  COVERED 


QUANTUM  MECHANICS  APPLIED  TO  DYE  LASERS 


17.  AUTHORO) 


Technical  Report 


6.  PERFORMING  ORG.  REPORT  NUMBER 


fa.  contract  or  grant  number^; 


io.  program  element,  project,  task 

AREA  &  WORK  UNIT  NUMBERS 

61102A  A3IB  JO 

i 

12.  REPORT  DATE 

September  1981 

13.  NUMBER  OF  PACES  ’ 

33 

15.  SECURITY  CLASS,  (of  th/e  report) 

UNCLASSIFIED 

"~15<».  DECLASSIFY  ATI  On7  DOWNGRADING 
SCHEDULE 


Carol  K.  Pearce 


9-  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

US  Army  Electronics  Research  &  Development  Command; 
Night  Vision  &  Electro-Optics  Laboratory; 

Fort  Belvoir,  VA  22060 

M.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

US  Army  Electronics  &  Development  Command  ;  Night  Vision 
&  Electro-Optics  Laboratory,  ATTN.  DELNV-L; 

Fort  Belvoir,  VA  22060 

M.  MONITORING  AGENCY  NAME  ft  ADDRESS (If  different  from  Controlling  Office) 


I  ie.  distribution  statement  (ot  thie  Report) 


Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (ot  the  abatract  entered  In  Block  20,  it  different  from  Report) 


Its.  supplementary  notes 


19-  KEY  WORDS  ('Continue  on  reveree  aide  If  nece emery  and  Identity  by  block  number ) 

Dye  Lasers  Emission 

Quantum  Mechanics 

Dyes 

Energy  Levels 
^^Njibsorption 

2<L~j^pTRACT  (Oothw  m  ravarem  a 7Se  tf  rtmeemeary  mod  Identify  by  block  number) 

The  application  of  quantum  mechanics  to  dye  laser  theory  is  reviewed,  based  upon 
recent  semiclassical  and  fully  quantum  mechanical  studies.  Both  semiclassical  and  fully 
quantum  mechanical  approaches  were  included.  Phenomena  investigated  in  the  dye  laser 
medium  included  spectral  hole  burning,  propagation  of  light  pulses,  and  light  absorption 
and  emission.  Some  studies  realistically  included  consideration  of  the  triplet  state.  Mix¬ 
tures  of  lasing  dyes,  or  of  lasing  and  saturable  absorber  dyes,  were  not  treated,  nor  were 
vapor  dye  lasers.  Theoretical  limitations  that  might  be  inherent  in  dye  lasers,  such  as 
|  gain  and  power  output,  are  other  important  areas  that  deserve  examination  in  future  studies. 


EDITION  OP  *  NOV  «5  IS  OBSOLETE 


_ UNCLASSIFIED _ \ 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  f*h«n  Date  Entered) 


CONTENTS 


Section 

I 

II 

III 

IV 

V 

VI 


I 

L 


Title  Page 

INTRODUCTION  AND  BACKGROUND  1 

APPLICATION  OF  THE  LAMB  LASER  THEORY  TO  6 

DYE  LASERS 

PROPAGATION  OF  LIGHT  PULSES  1 0 

SPECTRAL  CROSS-RELAXATION  AND  COLLISIONAL  1 2 

DEPHASING  EFFECTS 

A  QUANTUM-MECHANICAL  DYE  LASER  THEORY  1 6 

INVOLVING  EXACT  TREATMENT  OF  TRIPLET  STATES 

CONCLUSIONS  24 


Er- 

pj :  '  • 

l  v:  : 


PA 


if 


l 


Accession  For 

NT  IS  GRAS:  I  ~ 

STIC  TAP 

Unr  r»n  o  u:*i  c  o  -3  r~) 

Juc tif  iccr  o;: 


« 


lii 


QUANTUM  MECHANICS  APPLIED  TO  DYE  LASERS 


I.  INTRODUCTION  AND  BACKGROUND 

Successful  lasing  in  an  organic  dye  was  not  reported  until  1966  although  the  first 
successful  laser  was  developed  in  1961.  Use  of  dyes  as  saturable  absorber  laser  Q-switches 
preceded  the  laser  dye  success  and  continues  to  the  present.  NV&EOL  has  successfully 
applied  the  incorporation  of  a  dye  Q-switch  in  low-energy  laser  rangefinding  systems. 
High-energy  applications  were  explored.  A  concerted  effort  was  made  to  extend  the  long 
wavelength  applicability  of  the  dye  Q-switch  through  the  creation  of  new  dyes.  A  close 
relationship  exists  between  dye  lasers  and  dye  Q-switches.  One  of  the  first  instances  of 
lasing  in  a  dye  occurred  inadvertently  during  a  search  for  new  dye  Q-switches. 
Developments  in  the  more  widely  studied  phenomenon  of  dye  lasers  support  our  mission 
in  dye  Q-switches  and  have  provided  a  stimulus  for  this  report. 

Both  pulsed  and  continuous  wave  (CW)  dye  lasers  are  operational  in  configurations 
much  like  those  used  for  other  lasers  (except  for  modifications  required  to  offset  losses 
due  to  the  triplet  state).  Numerous  organic  dyes  have  been  found  to  lase;  the  cyanine, 
coumarin,  xanthene,  and  oxazine  classes  are  of  major  importance.1  The  tunability  in¬ 
herent  in  dye  lasers  is  probably  their  most  useful  asset.  Depending  upon  the  choice  of  dye, 
tunable  radiation  can  be  generated  from  the  near  ultraviolet  to  the  visible  and  near  in¬ 
frared  regions  of  light.  Linewidths  as  narrow  as  I0"4 *  nanometers  are  available,  as  are 
pulse  widths  of  0.01  nanosecond  duration.2  When  saturable  absorber  mode  locking  dyes 
are  incorporated  in  the  laser  cavity,  picosecond  pulsewidths  become  attainable.  Both  CW 
and  pulsed  (0.1  to  over  300  p/s)  operation  can  be  obtained  commercially,  with  peak 
power  up  to  2  x  106  watts  being  cited.2  Energies  between  several  microjoules  and  up  to 
over  10  joules  per  pulse  have  been  reported.3  The  potential  of  the  dye  laser  in  integrated 
optical  elements  has  been  discussed  with  possible  applications  in  communications  and 
data  transmission.3  The  precise  tunability  of  dye  lasers  mandate  their  use  in  spectroscopic 
studies  and  in  such  practical  applications  as  LIDAR.  LIDAR  (Light  Detection  and  Rang¬ 
ing)  enables  one  to  detect  and  measure  trace  atmospheric  pollutants  by  a  radar  type  of 
fluorescence  technique.3  Dye  lasers  have  been  investigated  for  optical  waming/detections 
systems  and  as  human  optical  countermeasure  sources.4  5  More  extensive  application  of 

*  K.  H.  Drexhage,  “Structure  and  Properties  of  Laser  Dyes,”  in  “Dye  Lasers,”  Chap.  4,  ed.  by  F.  P.  Schaefer, 
Springer-Verlag,  New  York  (1973). 

2  ' 

“Laser  Focus  Buyer's  Guide,”  pp  143-144,  Advanced  Technology  Publications,  Inc.,  Newton,  MA  (1979). 

3  tt 

T.  W.  Hansch,  “Applications  of  Dye  Lasers,”  in  “Dye  Lasers,”  Chap.  5,  cd.  by  F.  P.  Schaefer,  Springer-Verlag, 
New  York  (1973). 

4 

L.  Gampel,  S.  Guch,  and  H.  Wang,  “Dye  Laser  Module,”  U.S.  Army  Research  and  Development  Report,  ECOM 
76-1972-1  (Feb  77). 

^  D.  H.  Sliney  and  T.  L.  Lyon,  “Evaluation  of  a  Phase-R  Dye  Laser  in  Phase  II  of  an  ECOM  Countermeasure 
Source  Evaluation,”  AD-917  016L,  Defense  Documentation  Center,  Alexandria,  VA  (December  1973). 
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the  dye  laser  can  be  expected  in  all  of  these  areas,  particularly  as  the  state-of-the-art  ad¬ 
vances.  Some  of  the  improvements  that  should  develop  are:  dyes  and  host  materials 
designed  specifically  for  dye  laser  systems,  vapor  dye  lasers  operable  in  the  visible  and 
near  infrared  (and  pumped  by  electrical  discharge  or  electron  beam  sources,  and 
modifications  in  the  laser  module  itself. 

As  dye  laser  technology  requires  increasingly  unique  responses,  rather  than  merely  bor¬ 
rowing  from  the  existing  technologies  of  dyes,  photochemistry,  and  lasers,  guidance  will 
l>e  needed  from  theoretical  treatments.  Much  of  the  theoretical  work  done  on  dye  lasers  is 
based  upon  energy  balance  and  rate  studies.  Quantum  mechanical  treatments  of  dye  laser 
systems  are,  at  present,  much  less  numerous  than  for  other  laser  systems.  Among  the  few 
quantum  mechanical  dye  laser  studies  are  the  semiclassical  theories  of  Eberly,  Chin  and 
Zardecki,  and  Vahey  and  Yariv.6  7  8  A  fully  quantum-mechanical  study  was  reported 
recently  by  Schaefer  and  Willis.9  The  intent  of  this  report  is  to  discuss  dye  lasers  as  a  case 
of  light  interacting  with  matter,  the  case  being  treated  quantum  mechanically.  How 
quantum  mechanics  is  applied,  and  the  kinds  of  problems  solved,  will  be  of  major 
interest. 

Aspects  of  dye  lasers  covered  in  these  quantum  mechanical  treatments  to  be  discussed 
more  fully  below,  include: 

a.  Application  of  the  Lamb  laser  theory  to  dye  lasers.6 

b.  The  propagation  and  superposition  of  light  pulses  in  a  dye-solvent  medium. 

c.  The  absorption  of  light  in  organic  dyes  (both  in  dye  lasers  and  as  saturable  ab¬ 
sorbers)  as  it  is  affected  by  spectral  cross-relaxation  and  collisional  dephasing.8 

d.  Development  of  a  completely  quantum  mechanical  theory  of  dye  lasere,  including 
an  exact  calculation  for  the  dye  triplet  state.9 

Basic  to  all  of  these  studies  is  the  fact  that  organic  dye  molecules,  upon  absorption  of 
light  quanta,  undergo  an  optical  transition  from  a  ground  (singlet)  state  energy  level  to  a 
higher  excited  singlet  energy  level.  Decay  to  the  ground  state  results  in  fluorescent  emis¬ 
sion,  comparable  to  the  laser  emission  when  the  dye  is  in  a  laser  cavity.  Intersystem  cross¬ 
ing  (where  kST  is  the  intersystem  crossing  rate)  can  easily  occur,  however,  so  that  the  ex¬ 
cited  singlet  may  decay  instead  to  a  triplet  level  (Figure  1).  Thus  the  efficiency  of  the 
(singlet)  emission  process  is  reduced. 

6  J.  H.  Eberly.  “Toward  a  Theory  of  the  Dye  Laser,"  “Atomic  and  Molecular  Optics/*  Proc.  Rochester  Symp., 
ed.  by  J.  H.  Eberly,  Univ  of  Rochester  (March  1971);  also  US  Gov't  Publication  AD725218  (Defense  Docu¬ 
mentation  Ctr),  Alexandria,  VA. 

7  A.  L  Chin  and  A.  Zardecki.  Phyi.  Rev.  A,  13.  1528  (1976). 

8  D.  W.  Vahey  and  A.  Yariv.  Rhys.  Rev.  A,  10  (6)  1578  (1974). 

9  R.  B.  Schaefer  and  C.  R.  Willis.  Phys.  Rev.  A,  13  (5)  1874  (1976). 
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The  distribution  of  molecules  upon  interaction  with  light  will  be  among  a  minimum  of 
four  energy  levels,  so  that  the  organic  dye  laser  must  be  considered  as  a  four-level  system 
at  least,  in  contrast  to  the  two-level  laser  of  Lamb's  theory,6  Unlike  other  laser  media,  fur^ 
thermore,  organic  dyes  are  characterized  by  broad  absorption  and  emission  bands.  These 
broad  bands  afford  the  tunability  for  which  dye  lasers  are  known.  They  indicate  also  that 
many  vibrational/rotational  levels  are  populated  for  a  particular  electronic  transition. 

In  the  discussions  to  follow,  the  use  of  density  matrices  will  be  found  to  be  most  useful 
in  formulating  the  equations  of  motion  for  dye  molecules  moving  through  an  electric  field 
(the  field  due  to  the  electric  component  of  light  waves).10  11  The  density  matrix  is  one 
whose  elements  are  bilinear  combinations  of  the  probability  amplitudes,  CA  and  CB  (if 
^(r,  t)  is  a  known  wave  function):10 

t)  =  CA(#A(r)  +  C B(#B(r),  (1) 

where  \^(r)  is  the  stationary,  or  time-independent,  wave  function  for  molecules  A  or  B. 

Sargent,  Scully,  and  Lamb  (hereafter  SS&L)  have  discussed  the  above  “pure  case"  as  well 
as  the  more  generally  useful  “mixed  case"  when  the  state  vector  (in  Dirac  motion)  |  \^( t)> 
is  unknown.10  The  wave  function  ^(r,  t)  is  an  expansion  coefficient  of  the  state  vector.10 
Thus,  q  is  the  general  density  operator 


Q  =  ^  F  | 
j 


(2) 


P.  is  a  probability  factor  for  systems  which  have  a  state  vector  |  \£(t)>  and  where  two- 
level  state  vectors  are: 


and  so 


*,<»)>  =  Caj(t)|a>  +  Cbj(t)  |  b  > , 


|C  .  * 

1  aj 


K,!2/ 


(3) 


(4) 


^  J.  H.  Eberly,  “Toward  a  Theory  of  the  Dye  Laser,”  “Atomic  and  Molecular  Optics,”  Proc.  Rochester  Symp.. 
ed,  by  J.  H.  Eberly,  Univ  of  Rochester  (March  1971);  also  US  Gov’t  Publication  AD725218  (Defense  Docu¬ 
mentation  Ctr),  Alexandria,  VA. 

M,  Sargent  III,  M.  O.  Scully,  and  W.  E.  Lamb,  Jr.,  “Laser  Physics,”  Chaps  7,  14,  and  16,  Addison-Wesley,  New 
York  (1974). 

11  W.  E.  Lamb,  Jr.,  Phys.  Rev.,  134  A1429  (1964). 
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The  equation  of  motion,  utilizing  Schrodingeris  time-dependent  equation 

H>(r.t)  =  i1li*.(r,t)  (5‘ 

dt 

is: 

e  = -—  lH,ej.  (6» 

a 

SS&L  describe  the  use  of  the  density  matrix  equations  of  motion  in  various  aspects  of 
laser  quantum  mechanical  problems.10 

SS&L  showed  that  these  matrix  equations  of  motion  could  be  used  to  derive  equations 
for  the  Bloch  model.  In  the  Bloch  model,  a  two-level  system  of  magnetic  dipoles  of  spin 
±  Ms,  is  used  to  represent  another  two-level  system  (as  a  two-level  laser  system).  This  ap¬ 
proach  was  used  by  Vahey  and  Yariv®  for  the  singlet-triplet  dye  system.  The  Bloch  equa¬ 
tions  are  based  upon  the  equations  of  motion  through  a  magnetic  field  and  include  two 
relaxation  terms  of  equations.®  10  12 

Reservoir  theory  can  be  used  with  the  density  operator,  q .  in  laser  calculations.  The 
quantum  mechanical  treatment  of  dye  lasers  by  Schaefer  and  Willis9  combines  reservoir 
theory  with  density  matrices.  As  described  by  SS&L,  the  reservoir  approach  facilitates  the 
division  of  a  system  into  portions  of  primary  (A)  and  secondary  (B)  interest.  The 
“secondary”  part  becomes  the  reservoir,  which  can  be  eliminated  by  use  of  a  “reduced 
density  operator,”  (?AB(t).10 

The  “rotating  wave  approximation”  has  been  used  below.6  7  9  This  approximation 
allows  one  to  eliminate  terms  where  (a?  +  v)  is  in  the  denominator,  or  eMuJ  +  is 
present.6  10  Since  co,  the  “matter”  wave  frequency  in 

\p(T,  t)  =  Cal/-a(r)  e~iUa'  +  Cb^h(r)  (7) 


®  J.  H.  Eberly,  “Toward  a  Theory  of  the  Dye  Laser,”  “Atomic  and  Molecular  Optics,"  Proc.  Rochester  Symp.. 
ed.  by  J,  H.  Eberly,  Univ  of  Rochester  (March  1971);  also  US  Gov’t  Publication  AD725218  (Defense  Docu¬ 
mentation  Ctr),  Alexandria,  VA. 

7  S.  L.  Chin  and  A.  Zardecki,  Phys.  Rev.  A,  IS.  1528  (1976). 

8  D.  W.  Vahey  and  A.  Yariv,  Phys,  Rev.  A,  7 0  (6)  1578  (1974). 

9  R.  B.  Schaefer  andC.  R.  Willis,  Phys.  Rev.  A,  IS  (5)  1874  (1976). 

^  M.  Sargent  III,  M.  O.  Scully,  and  W.  E.  Lamb,  Jr.,  “Laser  Physics,”  Chaps  7,  14,  and  16,  AddisonAVesley.  New 
York  (1974). 

A.  Yariv,  “Quantum  Electronics,”  Chap.  8,  John  Wiley  and  Sons,  Inc.,  New  York  (1967). 


5 


and  the  "field”  frequency  of 


E(t)  oc  cos^t),  (8) 

are  equal  under  near  resonance  conditions,  and  w  >  >  1,  (od  4-  p)  is  relatively  large.  For 
eMu)  +  rapid  oscillation  occurs,  so  that  such  terms  average  zero.6 

In  each  of  the  discussions  below,  the  symbolism  as  used  by  the  authors,  respectively,  has 
been  maintained  in  the  section  discussing  their  work.  Kauzmann  has  discussed  the  fun¬ 
damentals  of  the  quantum  mechanics  of  light  interacting  with  matter  and  time-dependent 
processes.13 

II.  APPLICATION  OF  THE  LAMB  LASER  THEORY  TO  DYE  LASERS 

Eberly  sought  to  derive  a  dye  laser  theory  from  the  jointly  interacting  laser  electric  field 
and  the  polarization  density  of  the  dye.6  He  applied  the  Lamb  laser  theory  to  organic  dye 
lasers,  noting  two  important  modifications  would  be  necessary.11  The  dye  laser  must  be 
treated  as  a  four-level  system  at  least,  not  a  two-level  one,  and  Lamb’s  use  of  the  extreme 
Doppler  limit  would  be  inappropriate  since  dyes  are  conventionally  used  in  solution. 
(One  should  note,  however,  that  more  recently  lasing  has  been  achieved  with  organic  dye 
vapors.)14  15 

The  polarization  density  P(t,  z)  was  derived,  based  upon  the  Lamb  theory,  for  dye 
molecules  in  the  presence  of  the  perturbing  field  E(t,  z).  For  a  single  mode  of  oscillation 
of  an  electromagnetic  (light)  wave  in  the  laser  cavity,  given  by  U(z),  the  natural  frequency 
was  fl  =  v  4*  <}>  (where  <£(t)  is  the  function  describing  the  slow  variation  of  <f>  with  time, 
and  v  is  the  frequency  of  the  light  wave.) 

Further,  assumptions  were  made  that  discrete  modes  were  present  in  the  laser  cavity, 
that  the  vector  nature  of  the  field  could  be  ignored  (as  well  as  its  polarization),  and  that 
all  parameters  are  functions  of  time,  t,  and  cavity  axial  position,  z,  only.  For  the  perturb¬ 
ing  electric  field. 


E(t,  z)  =  U(z)  e  (t)  cos  (vt  +  <£(t)). 


®  J.  H.  Eberly,  ‘Toward  a  Theory  of  the  Dye  Laser,”  “Atomic  and  Molecular  Optics,”  Proc.  Rochester  Symp., 
ed.  by  J.  H.  Eberly,  Univ  of  Rochester  (March  1971);  also  US  Gov’t  Publication  AD725218  (Defense  Docu¬ 
mentation  Ctr),  Alexandria,  VA. 

11  W.  E.  Lamb.  Jr.,  Phys.  Rev.,  134  A 1429  (1964). 

W.  Kauzmann,  “Quantum  Chemistry,”  Part  V,  Academic  Press,  New  York  (1957). 

14  P.  W.  Smith,  P.  I\  Liao,  ct  at.,  Appl.  Phys.  Ltrs.,  25  144  (1974). 

^  B.  Steyer  and  !'.  P.  Schaefer,  Opt.  Comm.,  10  219  (1974). 


6 


Here  €  is  the  amplitude  and  is  assumed  to  vary  slowly  with  time.  Only  the  first  derivatives 
of  €(t)  and  <t>{t)  were  used  in  Eberly ’s  derivation.  Since  both  </>(t)  and  f(t)  have  been 
assumed  to  exhibit  little  change  with  time,  this  must  he  true  for  E(U  z)  also,  that  the 
treatment  is  semiclassical. 


The  time-dependent  projection  of  the  polarization  density  P(t,  z)  can  he  written  as  a 
function  composed  of  oscillations  in  and  out  of  phase  with  cost  (yt  +  of  E(t,  z): 


P(t)  = 


dz  U(z)  P(t,  z) 


cavity 


=  C(t)  cos^t  +  0)  +  S(t)  sin(yt  +  0). 


(9) 


Using  Maxwell's  wave  equation  modified  with  an  arbitrary  conductivity  constant  o  as  a 
cavity  loss  correction  factor,  coupled  differential  equations  for  e(t)  and  0(t)  were  obtained 
(where  17  is  an  arbitrary  constant  undefined  further  by  Eberly): 

(v  +  <£  -  Q)  e(t)  =  -fc  (vh)2) Qt)  (10) 

e(t)  +  Vi{olr]2)e{ t)  =  -^(Wt/2)  S(t).  (11) 

These  two  equations,  along  with  that  for  the  polarization  density,  P(t.  z),  are  those  equa¬ 
tions  subject  to  self-consistent  (field)  solution  in  EberJy’s  theory, 

Eberly  determined  the  polarization  density  according  to  Lamb's  laser  theory  but 
adapted  it  to  the  unique  characteristics  of  a  dye  system.  The  optical  transitions  involved 
in  a  lasing  dye  are  shown  in  Figure  1.  The  assumptions  basic  to  Eberly 's  treatment  are: 

a.  None  but  optical  transitions  contribute  to  the  polarization  density,  and 


n 

b.  P2(t,z)  =  £p*(tj,z) 


(12) 


where  n  =  3  in  the  present  rase  and  p.  is  the  partial  polarization  associated  with  the  i'th 
optical  transition. 

Addressing  oneself  to  one  of  the  transitions  involved  in  the  above  summation,  where 
n  =  (level  B  —  level  a),  and  specifying  XR(t)  =  rate-time'l-volume‘l  for  the  excitation  of 
molecules  into  level/state  “B."  P.(t)  becomes  PB(t)  or 


pBo)  = 


dto 


(13) 


7 


For  the  dye  molecules  undergoing  the  transition  B-a,  there  will  be  more  than  one  final 
state  possible,  as  noted  above,  so  that: 


=  £  (#iBa  8b  +  e£/ilB>.  (H) 

a 


where  ^Ba  and  /xaB  are  the  dipole  moments  corresponding  to  '.he  different  final  states, 
7tB(t,  tj  is  the  average  or  expectation  value  of  f*molecu[r  in  state  9  at  time  ^  and  Qti  is  the 
density’  matrix  of  that  same  molecule,6  J1 


From  the  equation  of  motion  for  this  density  matrix. 


(15) 


i  h  q  =  [H,  e)  +  JL  (70  +  GT)' 

7  being  the  diagonal  matrix:10 


(16) 


(17) 


where  7B  and  y  are  the  radiative  decay  constants  for  states  B  and  a,  and  following  the 
development  of  Eberly  and  Lamb,  and  invoking  the  rotating  wave  approximation,  the 
polarization  density  for  the  B— a  optical  transition  will  be:6  11 


(18) 


P„(t,  z)  - 


XB(t) 

^7 


L(~- 

•  V  h , 


l*Ba  **aB 


U(z)  e(t)e_i(vt  +  0) 


Ba 


v) 


+  complex  conjugate. 


J.  H.  Eberly,  “Toward  a  Theory  of  the  Dye  Laser,”  “Atomic  and  Molecular  Optics,”  Proc.  Rochester  Symp., 
ed.  by  J.  H,  Eberly,  Unjv  of  Rochester  (March  1971);  also  US  Gov’t  Publication  AD725218  (Defense  Docu¬ 
mentation  Ctr),  Alexandria,  VA. 

M.  Sargent  III,  M.  O.  Scully,  and  W.  E.  Lamb,  Jr.,  “Laser  Physics,”  Chaps  7,  14,  and  16,  Addison-Wcsley,  New 
York  (1974). 

1 1  W.  E.  Lamb,  Jr.,  Phys.  Rev.,  134  A 1429  (1 964). 


8 


Expressions  for  CH(t)  and  SB(t)  were  obtained  by  Eberly  from  the  P(t)  equation  (9): 


CB(t) 

al 

/  ^bn 

\  /3e(t) 

(19) 

V 

\7b, 

1 

V‘>  = 

-1 

(\\ 
l  A  J 

ae{\), 

(20) 

V  1 

where 


a  +  i  0 


MBaDa  ^Ba>d0JBa 
^Ba  +  i(tJBa  "  v) 


(21) 


here,  M  =  1,  since  restricted  to  single  mode,  and  Da(o>Ba)  =  energy  density  of  state  a  with 
energy  ^ia>Ba. 

Eberly  obtained  one  last  significant  result  from  the  theory  outlined  above,  for  the  zero 
order.  Using  the  above  expressions  for  CB(t)  and  SR(t),  the  photon  gain  from  the  B*—a 
transitions,  Gu  is: 

n  —  a 


G 


B 


<2e/e)B*.a  =  JL  SB(t)/e(t) 

V2 


MBa  Da  d  00 
[(CJ  -  l>)2  +  7Ba  ] 


(22) 


(23) 


This  expression  was  similar  to  one  Eberly  had  obtained  earlier  using  rate  equations  and 
applied  to  problems  assoeiated  with  the  threshold  for  organic  dye  lasers.4  14 

^  J.  H.  Eberly,  “Toward  a  Theory  of  the  Dye  Laser,”  “Atomic  and  Molecular  Optics,”  Proc.  Rochester  Symp., 
ed.  by  J.  H.  Eberly,  Univ  of  Rochester  (March  1971);  also  US  Gov’t  Publication  AD725218  (Defense  Docu¬ 
mentation  Ctr),  Alexandria,  VA. 

16  O.  G.  Peterson,  J.  P.  Webb,  W.  C.  McColgin,  and  J.  H.  Eberly,  J.  AppL  Phys.,42(5)  1917  (1971). 


III.  PROPAGATION  OF  LIGHT  PULSES 


The  semiclassical  self-consistent  approach  was  also  used  by  Chin  and  Zardecki  to 
describe  the  propagation  of  light  pulse  through  a  dye  medium  and  also  the  superposition 
of  short  light  pulses  in  that  medium.7  A  multi-level  model  for  a  dye  molecule  in  solution 
was  presented.  Triplet  states  are  ignored  since  the  duration  of  the  pulse  trains  is  too  short 
for  intersystem  crossing  to  occur.  Essentially,  the  “supermolecule"  (the  dye  molecule  in  its 
solvent  cage)  when  perturbed  by  the  absorption  of  a  light  pulse  undergoes  several  transi¬ 
tions  (A  — b,  b— B,  B-*a,  and  a— A  of  Figure  1)  between  equilibrium  (A  and  B  of 
Figure  1)  and  non-equilibrium  states  (b  and  a  of  Figure  1). 

Equations  of  motion  were  derived  using  a  population  matrix  formulation,  q(z,  t),  for 
the  interaction  of  the  radiation  field  and  the  dye  “supermolecules."  A  basic  assumption  is 
made  that: 


V’  <  <  T,,  <  <  7;'.  (24) 

where  Tp  is  the  pulse  duration,  yA‘l  is  the  atomic  memory  time,  and  7^  is  the  non- 
radiative  decay  time  from  b  to  B  in  Figure  1.  The  field  is  regarded  as  a  linearly  polarized 
running  wave  for  which  the  amplitude  and  phase  vary  slowly.  The  amplitude  envelope 
function,  £,  used  was: 


+  -L  =  -k£  -_L_  kSAB,  (25) 

3z  c  3t  2e0 

WL  mean  frequency  of  incident  pulse 

where  k  =  -  —  - *  K  represents  nonresonant 

c  velocity  of  light 

losses  in  the  medium,  and  €  is  the  vacuum  permittivity.  SAB  is  the  quadrature  component 
of  the  total  polarization,  PAR,  and  is  computed  using  the  rotating  wave  approximation. 
The  phase  function,  <£,  for  the  electric  field  is: 

<t>  =  co,  t  —  kz  +  0  (z,  t).  (26) 

As  in  Eberly’s  treatment,  a  perturbation  method  was  used,  the  perturbation  being 
/iE 

“ —  where  is  the  dipole  matrix  element  between  states  and  E  is  the  energy  of  the  transi  - 

tion.  Expressions  for  the  perturbed  and  unperturbed  Hamiltonians  were  given  to  be  used 
in  generating  the  equations  of  motion. 

7  S.  L.  Chin  and  A.  Zardecki.  Phys.  Rev.  A,  IS,  1 528  ( 1 976). 
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The  total  polarization  PAB  can  be  expressed  as  the  sum  of  the  partial  polarizations,  P^  : 


PA»  -  2  £  P, 

a  b 


(27) 


=  CAB  (z,  t)co8  <f>  +  SAB  (z,  t)sin  (f>. 


(28) 


Pab(Z'  *)  =  Aba  Pab  (Z-  *>  +  Mab  6ab  <Z>  *>• 


(29) 


The  population  maxtrix  used  to  express  the  equations  of  motion  for  the 
“supermolecule”  in  any  of  the  states,  a,  in  levels  A  or  B  was: 


e(z,  o 


E 


dtnX„<Z-  U  C(«*  Z'  ‘n- 


(30) 


Here,  p(a,  z,  t  t)  is  the  density  matrix  for  a  molecule  (at  place  z,  time  t#)  in  state  a,  and  X 
(z,  tj  is  the  molecular  concentration  per  unit  time  in  state  a ,  Four  equations  of  motion. 


-i  ^  *  ,  ,  . 

Q  =  -  (H,  Q),  were  given,  gaa<,  Qhh  ^  gha,  and  £ah.  Using  zeroth  order  elements  of  the 

it 

population  matrix.  Chin  and  Zardecki  define  the  population  difference  in  a  and  b,  : 


Nff  =  9S 


(31) 


If  the  photon  flux  is  defined  as  I  =  ’  and  In  is  the  input  pulse  at 

z  =  0,  then  the  flux  modified  by  multiple  transitions  is: 


1(2,  t) 


_ Ip  (t  —  z/c)  e(a-  -  2k)? _ 

1  [a2/(a,  -  2k)]  I0(t  7./c)  {l  e l,a‘ 


(32) 


Resonant  absorption  to  the  first  and  second  order  in  the  photon  flux  is  given  bv  the  coeffi¬ 
cients  at  and  ar 


II 


Chin  and  Zardeeki  applied  their  theory  to  the  spectral  hole  burning  phenomena,  ob¬ 
taining  results  consistent  with  those  of  Mouro17  for  dye  molecules  in  a  solvent  cage.  They 
also  addressed  the  case  of  two  adjacent  pulses  being  superposed  in  a  dye  medium.  The 
first  pulse  obeys  the  condition  of  equation  (24).  The  second  pulse  is  similar  to  the  first 
and  occurs  after  a  time  delay,  T.  This  delay  is  much  shorter  than  the  fluorescence  time  of 
the  dye  molecule  tit  level  B  and  is  much  shorter  than  the  solvent-solute  relaxation  time.  As 
in  the  propagation  study,  equations  of  motion  were  derived  for  the  time  evolution  of  the 
population  matrix  elements,  and  expressions  for  the  photon  flux  were  found.  The  pulse 
superposition  effect  is  defined  by  the  authors  as  the  storage  of  the  energy  of  a  first  pulse 
as  a  population  inversion,  and  stimulated  emission  occurring  subsequently  upon  irradia¬ 
tion  with  a  second  pulse.  Their  derivation  involved  separating  the  absorption  and 
amplification  processes,  so  that  in  the  direct-transition  approximation  the  photon  flux 
becomes: 


31  +  I  31 
3z  c  3t 


eBB<T )Zom-2k-  1+^1  I 


I  (33) 


Qualitative  experimental  support  for  the  superposition  theory  was  presented  using 
oxazine  perchlorate  in  glycerine. 

IV.  SPECTRAL  CROSS-RELAXATION  ANI)  COLLISION  AL  DEPHASING  EFFECTS 


That  the  absorption  of  light  by  organic  dye  molecules  is  far  from  a  simple  process  was 
shown  in  the  scmielassical  theory  of  Vahey  and  Yariv.8  They  treated  dye  solutions  in¬ 
teracting  with  light  at  high-incident  flux,  as  an  ensemble  of  optical  dipoles.  Absorption  of 
light  in  this  case  mav  exhibit  inhomogeneous  broadening  or  absorption  occurring  for  the 
ensemble  over  a  wider  frequency  band  than  would  be  anticipated  for  a  lone  optical 
dipole.8  10  Spectral  cross-relaxation  results  when  the  configuration  of  the  solvent  cage  or 
environment  surrounding  the  optical  dipole-dye  molecule  is  altered  in  some  way.  Another 
effect  which  mav  occur  during  the  interaction  of  light  and  an  optical  dipole-dye 
molecule,  is  that  of  a  random  shift  in  the  phase  of  the  radiating  dipole  due  to  collisions  in 
the  dye  system.8  10 

8  D.  W.  Vahtfy  and  A.  Yariv,  Phys.  Rev.  A,  10  (6)  1578  (1974). 

^  M.  Sargent  III,  M.  O.  Scully,  and  W.  F.  Lamb,  Jr.,  “Laser  Physics,”  Chaps  7,  14,  and  16,  Addison-Wesley,  New 
York  (1974). 

17  Ci.  Mouro.  IM  F  J.  Quantum  Electron..  QK  l /  (I)  l  <t975). 
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Vahey  and  Yariv  point  out  that  earlier  treatments  of  the  interaction  between  light  and 
matter  have  included  neither  spectral  cross-relaxation  effects  nor  collisional  dephasing. 
They  presented  a  semiclassical  study  which  did  consider  these  effects.8  Application  of 
their  theory  to  experimental  data  for  certain  well-known  cyanine  saturable  absorber  dyes 
accounted  for  observed  ‘*hole-bumingM  in  the  absorption  spectrum  due  to  cross-relaxation 
(see  Figure  3  of  Mouro17).  They  determined  that  if: 

Tj  =  time  constant  for  deactivation  or  decay  or  excited  molecules 

T2  =  time  constant  for  collisional  dephasing 

T3  =  time  constant  for  spectral  cross-relaxation  (and  if  TA  =  time  constant  for  absorp¬ 
tion,  and  —  oc  rate  of  optical  process),  then  distortion  in  the  absorption  spectrum 

will  occur  when  —  >  — ,  and  coherent  interaction  effects  will  be  significant 
Ta  T3 

when  —  ~  1—  Vahey  and  Yariv  are  concerned  here  with  systems  in  which 

Ta  T2 

1 

T2  <  “  Tr  Rate  equations  were  derived  by  Vahey  and  Yariv,  based  upon  concepts 

evolving  from  likening  the  dye  absorbers  to  an  ensemble  of  liquid-free  spins  interacting 
with  radiation.  The  Bloch  equations  were  applied  to  such  systems,  affording  a  vector 
model  of  the  density  matrix  equations  of  motion.8  10  12  The  model  used  was  that  of  a 
single  dye  molecule  (state  vectorin',  v',  r'))  which  underwent  only  one  collision,  a  strong 
collision,  at  time  to.  The  undamped  Bloch  equation. 


dr/dt  =  w  x  r' 


(34) 


was  applied  to  the  molecule  with  the  conditions  that: 

r'  (t  =  to)  =  (0,0, ±  1);  ir^(t)  =  (u\  v\  r') 
uT  =  (cj,(x,  t),  0,  a>o  —  co  +  6(t)) 

(a>  —  u)  )is  the  detuning  frequency, 

o)}  =  a  frequency  proportional  to  the  field  amplitude, 

6(t)  =  the  temporal  detuning  generated  by  cross-relaxation. 

8  D.  W.  Vahey  and  A.  Yariv,  Phys.  Rev.  A,  10  (6)  1578  (1974). 

M.  Sargent  III,  M.  O.  Scully,  and  W.  E.  Lamb,  Jr.,  ‘‘Laser  Physics,”  Chaps  7,  14,  and  16,  AddisonWcsIcy.  New 
York  (1974). 

^  A.  Yariv,  “Quantum  Electronics,”  Chap.  8,  John  Wiley  and  Sons,  Inc.,  New  York  (1967). 

^  G.  Mouro,  IEEE  J.  Quantum  Electron.,  QK-ll  (1)  1  (1975). 
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A  solution  of  the  Bloch  equation  is  v'(o)r  O)o  —to,  t-to,  dt  1).  Since  the  undamped  equa¬ 
tion  was  used,  Vahey  and  Yariv  note  that  u',  v',  and  r'  are  not  ensemble  averages,  but 
are  averages  determined  by  a  summation  of  microscopic  state  vectors  produced  since 
times  tu  by  random  collisions  at  interval  T2. 

^Using  this  model  to  address  the  question  of  optical  absorption  in  dye  molecules,  where 
p(x,  t,  ojJ  is  the  induced  dipole  moment  (due  to  a  polarized  field,  E(x,  t))  and  pu  is  the 
dipole  matrix  element  due  to  the  initial  and  final  states,  the  authors  obtain: 

— ♦  ■  # 

p(x,  t,  U>(>)  =  pj— u(x,  t,  w,)eos(w  t-kx)  +  v(x,  t,  oj(i)sin(u>t-kx)]  (35) 


since 


E(x,  t)  =  Pz  £  (x,  t)co6(oj  t-kx). 
The  out-of-phase  component  of  polarization,  v(x,  t,  u>J,  is: 


(36) 


v(x,  t,  toQ)  = 


/t  dt 

— ^  e-(t-to)/T2  (nb(x,  to+,  co0)  v'(W|,  w,  t  -  tfl 

T2 


1) 


(37) 


+  na(x’  to+’  W  v'(wl .  «0  “  CJ»t  ~  t0  +  1)1 . 


This  is  the  instantaneous  polarization.  The  probability  that  the  dye  ground  state  re¬ 
mained  occupied  after  collision  occurred  (at  time  to)  was  given  by  njx,  t  + ,  w  ).  The 
probability  that  dye  molecules  would  be  found  in  an  excited  state  at  time  (t  -I- )  was 
n  (x,  to  + ,  cjJ;  to  eliminate  systems  showing  stimulated  emissions,  na  =  0.  Development 
of  an  expression  for  n,(x,  to  +  ,  u>o)  leads  to  expressions  for  the  induced  polarization: 


V(x,  t,W0)  = 


_ C0i(x,  t)  T2 _ 

1  +  c0|  (x,  t)  T2  +  (to0  -  oj)2  T| 


|nb(x,t,«0)[l  -(T,/T3)l  (38) 


d  Wo  8(^0  -  w0)mb(x,  t,  a/ )(T2/T3)  +  nc(x. 


t)(T2/T,)| 


In  phase  «(x,  t,  a)J: 


U(X.  I.  Uj  =  -  (0>(1  -  «(  T2v(x,  t,  wj. 


(39) 
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To  relate  the  induced  polarization  and  light  absorption,  the  electromagnetic  wave  equa¬ 
tion  was  separated  into  two  equations,  one  for  phase  and  one  for  amplitude.  The  wave 
amplitude,  £(x,  t)  is 


—  +  —  —  =  ~'/2No - L  f  dw0g(w0  -  wo)  v(x,  t,  w0).  (40) 

3x  c  3t  €  n  c 

o  o 

Here,  No  is  the  absorber  concentration,  while  the  refractive  index  of  the  solvent  is  given 
by  nw.  One  of  the  basic  consequences  obtained  by  Vahey  and  Yariv,  the  photon-transport 
equation,  evolved  from  the  wave  amplitude  equation.  The  transport  equation  is: 

f—  +— —  F(x,  t)  =  -  NoahF(x,  t)<  n^x,  t,  w0)  £  (F,  wo  -  ojJ>  av„  (41) 

\9x  c  3t  / 

where  the  flux  density  is: 

F(x,  t)  s  V4  eo  cno  £2(x,  t)/hu>,  (42) 

the  homogeneous  cross  section  for  interaction  of  molecules  in  exact  resonance  with  the 
wave  is 


°h  =  ^Po  Tj/e0cn01t. 

The  effective  ground-state  population  is 

nff(x,UW(,)  =  n^x.t.ojJM^/T,)) 
+  <  \(x,UWo)  >av  (T2/T3) 

and  £  (F,  w  —  «)  is  a  modified  Lorentzian  function. 


(43) 


(44) 


The  other  significant  theoretiral  result  obtained  was  the  population  rate  equation: 

—  nb(x,  t,  w0)  =  -  on  F(x,  t)  neff  (x,  t,  w0)  JC(F,  u>o  -  w)  - 
9t 


{lnb(x,  t.  wD)  -  <  nb(x,  t,  w0)>  av)/Tj }  +  nc(x,0T, 


(45) 
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Vahey  and  Yariv  applied  their  theory  to  two  classes  of  absorbers,  according  to: 

(1)  Steady-state  limit,  where  the  pulse  duration  is  long  with  respect  to  T(  and  T.,. 

(2)  Short-pulse  limit,  where  the  pulse  duration  is  very  short  compared  to  Tj  and  T.,. 

For  the  steady-state  case,  they  found  that  the  experimental  parameter  in  pant-pulse 
studies  with  organic  dyes  is  the  steady-state  transmission: 


F(x>  =  Q2  (x)  -  1 
F(0)  Q2(0)  —  1 


(46) 


where 


Q  -  M  +  Fah(T2  +  T,)f  (47) 

and  F(0)  is  a  known  input  flux  deifr^.  In  trie  .'hod-pulse  case,  the  theory  w  as  extended  to 
such  problems  occurring  in  real-w* 'H  '<  absorbing  systems  as: 

(a)  Ex  i ted  state  absorption. 

(b)  The  molecular  anisotiopv  dependence. 

(c)  The  multiplicity  of  interacting  excited  states. 

V.  A  QUANTUM-MECHANICAL  DYE  LASER  THEORY  INVOLVING 
EXACT  TREATMENT  OF  TRIPLET  STATES 

Schaefer  and  Willis  were  the  first  to  present  a  dve  laser  theory  which  w  as  not  only  fully 
quantum-mechanical,  including  radiation  fluctuations,  but  also  treated  the  triple  states 
exactly.9  As  noted  above,  dye  lasers  are  unique  among  lasers  because  both  singlet  ami 
triplet  levels  can  be  populated  during  lasing. 

The  model  upon  which  this  quantum-mechanical  treatment  was  based,  consisted  of  two 
main  parts— matter  and  radiation.  Reservoir  theory  was  invoked,  creating  a  matter  or 
molecule  reservoir  (singlet  and  triplet  level  decay)  and  a  radiation  reservoir  (accounting 
for  radiative  losses  in  the  laser  cavity).10 

9  R.  B.  Schaefer  and  C.  R.  WiUij,  Phys.  Rev.  A.  /J (5)  1874  (1976). 

M.  Sargent  III.  M.  O.  Scully,  and  W.  K.  Lamb,  Jr.,  "Laser  Physics,"  Chaps  7.  14,  and  16.  Addison-Weslc).  New 

York  (1974) 
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The  Hamiltonian,  H,  for  this  model  became; 


H  =  IHK  +  UM  +  VK_J  +  |Hr_  +  VK_rJ  +  |Hm_  +  VM_  J.  (48) 

where  V  represents  the  intersection  Hamiltonian  potential.  The  terms  are  defined  by 
Schaefer  and  Willis  as; 


Hr  =  radiation  Hamiltonian 

Hm  =  molecule  Hamiltonian 

/n 

HR  ^  =  radiation-reservoir  Hamiltonian 

/N 

H„  =  molecule-reservoir  Hamiltonian 

VI— r*> 

VH  >I  =  radiation-molecule  interaction 
V„  =  radiation-radiation  reservoir  interaction 

Vvl  =  molecule-molecule  reservoir  interaction. 

>1— rt> 


The  reservoir  terms  were  eliminated  in  Schaefer  and  W  illis  because  the  only  reservoir  at¬ 
tributes  of  interest  were  those  of  pumping  and  decay  and  not  specific  details  about  the 
reservoirs. 


The  radiation  field  Hamiltonian  was:9 


HX  =  5  £  Ti  J2k  +  '/*) 


K,a  k,a 


(4»» 


where  -►  represent*  summation  over  all  cavity  modes  (in  ihe  present  ease  k  =  1. 


since  only  a  single-mode  was  treated),  and  ^  represents  summation  overall  a 
polarizations.  The  frequency  associated  with  k  is  The  “a  and  a*  M  operators  an' 
annihilation  ami  creation  photon  operators,  respectively.10 


9  R.  B.  Schaefer  and  C.  R.  Willis,  Phys.  Rev.  A,  U  (5i  1874  (1976) 

^  M.  Sargent  III,  M.  O.  Scully,  and  W.  K.  Lamb.  Jr.,  “Laser  Physics,"  Chaps  7.  14.  and  16.  Addison-Wesle\ .  Now 
York  (1974). 


Turning  to  the  Hamiltonian  for  matter,  HM,  the  singlet-triplet  states  associated  with  dye 
molecule  energy  levels  require  that: 

/\  /\ 

H„  =  HJ+  Ht  (50) 

where  _ 

ft.,  =  ?  ?  (so,i  VI  ei} 

'h't  =E  E<t£j  to  j  Cj) 

o  j 


and  o 


ei&ej 

Sl.i  &  Cj 


S„  &  t  . 
0,1  a,j 


The  relations 


=  total  dye-molecules  active  in  laser  system 

=  singlet  states 

=  triplet  states 

=  energies  corresponding  to  i  or  j  levels,  respectively 
=  creation  operators  for  singlets,  triplets,  respectively 
=  annihilation  operators  for  singlets  and  triplets,  respectively, 

hip  between  the  molecular  annihilation  and  creation  operators  were  given: 


1st 

0,1 

,  s  .  =  6] 

’  OJ  1JJ 

-  V 

(51) 

st 

a,r 

O  — F 

II 

Ivi-Vjl  * 

■•;.,>  • 

t«.j-*.jl =0 

(52) 

[sf 

0,1 

K 

IVrVi1  =  0 

(53) 

That  the  molecule  or  matter  terms  must  be  treated  as  singlets  and  triplets  was  shown  also 
in  the  radiation-matter  interaction  Hamiltonians: 


V„_>,  =  VK_S  +  VH_r  (54) 

From  llu*  matter-radiation  interaction  (dipole  approximation).  Schaefer  and  L  illis  found: 


(55) 


8 


(55) 


VR-S  =  *  53 


a 


53  L  so,j  4  +  so,i  V- 

i,j  k 


where  e.  >  e 

T,,,,!,  =  (4  ir  e*(Cj  -  Cj)2  hQk  Vf  x  (T •  x>  T  (x); 


VR-T  ~  *  53  Z)  13  y olmk^o.l  to,m  ak  +  ll.m  lo,l 
°  I*  m  k 


where  e,  >  e  , 

I  in 


(56) 


TW  =  (4  T  e3(fm  -  fi)a  h  vr  *  (X,J  r  ( X ). 


m  =  electron  mass 


X|m  =  dipole  matrix  element 


e  =  electronic  charge 


V  =  volume  of  laser  cavity 


r  (x)  =  eigenfunction  for  the  spacial  cavity 


£  for  j  >  i 
ij 

S  for  m  >  1 

1,  m 


I]  =  1  for  single-mode  operation. 

The  rotating  wave  approximation  was  applied  to  the  interaction  Hamiltonians  by  the 
authors.  (Thus  transitions  were  eliminated  where  photon  absorption  might  produce  a 
lower  energy  level  and  photon  emission  might  yield  a  higher  energy  level).  Thus,  from  the 
commutation  relations  cited,  s*  s  (or  t*  t  )  terms  correspond  to  the  creation 
(emission)  of  a  photon  and  decrease  in  molecular  energy  level.  Similarly,  s*  s  a(or  t*  t  a) 
terms  indicate  photon  annihilation  (absorption)  and  an  increase  in  molecular  energy 
levels. 
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fl 


.0 

Having  defined  the  system  Hamiltonian.  Schaefer  and  Lillis  utilized  the  Born- 
Markoff18  approximation  for  the  Liouville  equation. 


d'p'(t) 
d  t 


=  —  [H,^  (t)|,  p  (t)  being  the  density 

dr 


operator  of  the  system 


(57) 


To  demonstrate  the  validity  of  such  an  approach,  they  calculated  that  experimental  data 
for  Rhodamine  6G,  one  of  the^est-studied  laser  dyes,  did  indeed  satisfy  the  criterion  for 
validity  that  its  (continuous-wave  lasing)  interaction  time  is  minimal  compared  to  the 
relaxation  time,  or  (t.^/t^)  <  <  1. 

Thus, 


and 


U.025  </-S  N  \  <0.04 

V  S.i  / 


0.01  <^  TN^j  <0.03 


(58) 


(59) 


where 

N  =  number  of  lasing  organic  molecules 
S  =  stimulated  emission  rate  of  singlets 
T  =  stimulated  absorption  rate  for  triplets 
*s  ,  =  lower  singlet  decay  rate 
VT  u  =  uPPeT  triplet  decay  rate. 

Clearly,  both  the  singlet  and  triplet  contributions  satisfy  the  (t^/t^)  <  <  1  criterion. 

18  SS&L,  in  discussing  the  Markoff  approximation,  note  that  the  average  value  of  a  random  shift  in  an  energy 
difference  is  zero,  so  that  if  the  variations  of  this  shift  are  rapid  compared  to  other  decay  times  (relaxation 
,  (average  shift) 

times)  in  the  system,  then,  ■  ~  0. 

decay  time 
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'X  /X 

If  H  =  H  +  X  V  and  X  is  dimensionless,  the  above  approximation  was  given  as: 


+  _L  [  Ho,/p(t)|  =  -X  2  f  dTlV,(V(T)/fc], 
dt  h  J 


(60) 


where  V(T)  =  e,l,°T/h  V  e  lH°T/h. 


(61) 


This  equation  was  shown  by  Willis  to  describe  laser  photon  statistics  if  one  does  not 
separate  matter-radiation  interactions  into  factorable  matrices.19 

Equations  of  motion  were  derived  next.  As  derived  bv  Schaeffer  and  W  illis,  using  the 
Born  approximation  master  equation  where  the  diagonal  matrix  is  coupled  only  to  the 
diagonal  elements,  two  matrices  result: 


Radiation  Density  Matrix 
R  =  <  n I  R I  n> 

n  1  * 

Diagonal  Single-particle  Radiation  Density  Matrix 

Dni(t)  =  <  n |  <  j  |  D(l,  t)  |  j>  |  n>  ; 


(62) 


(63) 


here  n  is  the  photon  number,  and  j  refers  to  the  state  of  the  jth  molecule.  Two  sets  of  equa¬ 
tions  of  motion  were  formulated  from  these  matrices.  One  set,  the  equation  of  motion  for 
the  radiation  density  operator^t)  is: 


d  R  (t)  ©  © 

_ 2--  =  -<■+  i)Rn+,l 

dt 


©  © 

S  N  [(n  +  1)  (DS'U  -  D*:6! )  +  n(Dns’e  -  D*-“ ,  )1 


(64) 


T  N  [(n  +  1)  (Dj,u  -  Dj;Sj )  + 


19 


C.  R.  Willis,  Phys,  Rev.,  165  420  (1968). 


21 


The  terms  labeled  (^P)  and  (jP)  give  the  radiation  decay  associated  with  the  radiation 
reservoir,  (IT)  is  emission  (stimulated  and  spontaneous)  in  the  singlet  state,  (^P)  is 
stimulated  absorption  for  the  ground  singlet  state,  and  (^5^  and  are  the 

terms  corresponding  to  and  (^4^  respectively,  for  the  triplet  state.  The  other 


/\ 


set  for  the  singlet-molecule-radiation  density  operator  D(f,  t)  was  composed  of  four  equa¬ 
tions  to  satisfy  the  possibility  that  the  molecule  could  be  in  a  single  (S)  or  triplet  (T)  state, 
and  in  an  upper  (u)  or  lower  (£)  level  of  state: 


dDsn’u(t) 


dt 


^S.U  +  Ksr>  DnU  +  R„(t)  *e  -  S(n  +  1 )  (Ds/  -  Ds/, ),  (65) 


dDs/  (t) 
dt 


=  -'s..  ~  sn  «*»*  -DS.“l)(n  +  l) 


(66) 


dDT-u  (t) 


dt 


=  -^T  u  DT  “  -  T(DT.u  _  DT/,)(n+  1) 


(67) 


dDT/  (t) 


dt 


=  Kst  D».«  +  ^T,u  DT/_ ,  -  ,  -  n  1  (DT/  -  DT/_  , (68) 


These  two  sets  of  equations  are  the  fundamental  ones  for  the  Schaefer-Willis  dye  laser 
theory.  Several  solutions  were  obtained,  incorporating  the  parameters  a  and  (3: 

a  =  a  triplet  parameter: 

the  value  of  a  increases  with  triplet  concentration. 

a  m  (T  Kst/S  pT(). 
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0,  where  0  <  0  <  1,  and 


where 

T  =  triplet  stimulated  absorption  rate 
S  =  singlet  stimulated  emission  rate 
Kst  =  intersystem  crossing  rate 
0  =  fraction  return  of  excited  triplets,  0  <  0  <  1 
=  molecular  excitation  rate  due  to  external  pump 
n  =  number  of  photons 
VR  =  radiation  decay  rate. 

The  dye  laser  exact  solution  for  the  steady  state  was  obtained: 


^n+l=  ^p=0  G(P) 
where  G(p)  = 


St?Su  +  Tv 

P  P 


T.u 


+  Tt»I+ j  +^/NJ 


p  is  an  integer. 


I]  ^  is  the  product  of  terms  for  0  <  p  <  n. 


and  T?sn,u  is  obtained  from  D^,f  (t)  =  t}*'1  (t)  Rn  ( t >  when 
dDsnu  (t) 


dt 


=  Q  in  equations  (65)  —  (68). 
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Simplified  solutions  of  this  equation  were  obtained  by  (1)  eliminating  the  triplets,  giv¬ 
ing  the  "singlet-state  ease."  and  (2)  ineluding  the  triplets,  but  neglecting  terms  for  triplet- 
state  emission  and  singlet-state  absorption  to  obtain  the  "approximate  solution/'  Finally, 
graphical  solutions  of  the  exact  equation  were  obtained  by  computer  methods. 

Results  of  the  vari<  us  solutions  were  expressed  as  the  extrema  of  R|(,  the  threshold  value 
of  pumping,  and  variations  in  Rn  when  pumping  values  are  above  threshold.  One  in¬ 
teresting  result  of  this  theory  is  that,  instead  of  one  (maximum)  extremum  for  the  radia¬ 
tion  distribution.  R  ,  as  is  found  for  a  gas  laser,  for  the  dye  laser  it  is  possible  to  expect 
both  a  maximum  and  a  minimum  for  Rn,  depending  upon  other  laser  parameters. 

VI.  CONCLUSIONS 

The  basic  quantum  mechanical  studies  of  dye  lasers  have  been  reviewed.  Despite  the 
complexity  of  dye  laser  systems,  compared  to  other  laser  systems,  successful  results  have 
been  achieved.  An  acceptable  understanding  exists  of  dye  lasers  and  their  potential 
prospects  in  the  areas  of  spectral  hole  burning,  gain,  propagation  of  light  pulses  in  the 
dye  medium  (and  the  superposition  of  two  light  pulses),  light  absorption  and  emission 
processes  for  states  of  different  multiplicities,  and  pumping  thresholds.  Although  much 
remains  to  be  done,  these  studies  demonstrate  the  unique  contribution  quantum 
mechanics  can  make  upon  an  examination  of  a  time-dependent  interaction  of  matter  with 
the  electrical  field  component  of  light  in  a  dye  medium. 
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